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Abstract 

The new kind of events with a rapidity gap between two high-Er jets, observed in 
high energy pp collisions at the Tevatron, is found to be well described by the exchange 
of a colour singlet gluon system in the BFKL framework. This requires going beyond 
the conventional asymptotic Mueller- Tang approximation, which results in qualitatively 
different features of the basic parton-parton scattering amplitude. Non-leading corrections 
to the BFKL equation are included by incorporation of the consistency constraint and 
the running QCD coupling. Hadronisation and other non-perturbative QCD effects are 
treated through a complete Monte Carlo simulation, providing a gap survival probability 
that varies event-by-event, facilitating comparison with experimental results. 



1 Introduction 



The phenomenon of rapidity gaps, i.e. a region in rapidity or polar angle with no final state 
particles, has received much interest in recent years. This goes back to conventional diffractive 
scattering where a beam proton is quasi-elastically scattered, keeping most of its original mo- 
mentum (xp > 0.9), and emerges separated by a rapidity gap from the remaining final state 
particles. The introduction of a hard scale in the process |TJ has opened new possibilities to 
investigate an underlying parton process using perturbative QCD (pQCD). In spite of the hard 
scale, the rapidity gap is here always produced through soft, non-perturbative QCD dynamics 
as described in the pomeron model |I[ or the soft colour interaction model ||. This is given by 
the fact that the forward proton (or small mass system) is not affected by any large momentum 
transfer. 

A new class of rapidity gap events has recently been discovered in pp collisions at the 
Fermilab Tevatron ||, [|. The gap is then central in rapidity and has a high-E'r jet on each 
side corresponding to a large momentum transfer t ~ —E\^ et across the gap. The models for 
soft gap formation do not give a satisfactory description of this new phenomenon ||, which 
calls for an explanation in terms of a hard colour singlet (HCS) exchange in terms of pQCD. 
Mueller and Tang proposed the use of the BFKL equation |7| to describe the exchange of 
a colour singlet gluon ladder, as illustrated in Fig. |lja, giving a rapidity gap between the two 
hard-scattered partons. Their asymptotic solution does not, however, give a good description 
of the experimental data || f|. 

In this paper we provide a more complete, non-asymptotic solution to the BFKL equation 
for this process including also formally sub-leading effects. By implementing it in the Pythia 
Monte Carlo || we take into account the soft QCD processes that are necessary to produce the 
hadronic final state. The soft processes can destroy the parton level gap and result in a gap 
survival probability significantly smaller than unity. Our complete model gives a good descrip- 
tion of the observed 'jet-gap-jet' events. This provides an important test of the BFKL equation 
which describes novel QCD dynamics, beyond the conventional DGLAP parton evolution ||. 



2 The BFKL calculation 

The leading logarithmic (LL) BFKL equation [[j]] was introduced to describe the total cross- 
section for scattering of two small objects by summing, in the leading log(l/x) approximation, 
gluonic ladder diagrams with reggeised gluons along the ladder. The BFKL equation naturally 
generates the steep rise at small x of the gluon density in the proton and is of importance 
for the present exploration of QCD dynamics at HERA as well as for hadron colliders. The 
LL BFKL equation gives, however, a too steep energy dependence of the scattering amplitude 
and suffers from consistency problems due to the diffusion of gluon momenta into the infrared 
domain. Substantial progress in how to cure these problems has been achieved in the last 
few years. The next-to-leading logarithmic (NLL) corrections to the BFKL kernel have been 
computed [10] and, after a proper treatment of the collinear divergencies flllf , result in a 



lower intercept of the hard pomeron (about 1.2-1.3) which agrees with the experimental value. 
The infrared problems seem to be greatly reduced by unitarisation corrections to the BFKL 
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Figure 1: (a) The jet-gap-jet event topology in pp scattering resulting from hard colour singlet 
gluon ladder exchange between two partons as described by the BFKL equation. (b,c) The 
impact factors for the pomeron coupling to the same parton and to two different partons. 



pomeron [|T^, which suppress the evolution for low momenta by adding a non-linear term 
in the BFKL equation. The resulting phenomenon of saturation is confirmed by the success 
of the saturation model []14[. Thus, it is important to include these theoretical improvements 



when making critical tests of the BFKL equation against experimental data. 

The BFKL formalism is also capable of taking into account non-forward hard colour singlet 
exchange, which plays a crucial role in high-t elastic vector meson production [TjJ and in 
events with gaps between jets || The latter phenomenon is here studied in a theoretically 
improved BFKL treatment of the complex gluonic object constituting the hard pomeron. For 
asymptotically large rapidity separation y between the jets, it has been shown || [L7| that the 
LL BFKL pomeron may be regarded as coupling to individual partons resolved at the scale \t\. 
The origin of this factorisation property which should hold beyond the LL approximation, is 
motivated as follows. 

When the rapidity y is small the BFKL ladder (fig. [l|a) reduces to two simple gluons with 
transverse momenta ki^, with (k\ + A^) 2 — —t. In the two gluon approximation there appear 
contributions from scatterings on one parton (fig. |l]b) and two partons (fig. |l]c) in each proton. 
The contributions of individual diagrams involving the subdiagrams (figs. [l]bc) are infrared 
divergent when — > . However, in the sum of all the contributions these divergencies 
cancel due to colour coherence and the colour neutrality of the proton and the final state. 
When both the exchanged gluons couple to the same pair of partons (as in fig. |l]a), the loop 
integration gives a large logarithm of the momentum transfer while in the other case (when 
subdiagram fig. |l]c is relevant) the scattering amplitude is not logarithmically enhanced. For 
increasing y, when BFKL evolution becomes important, the contribution of gluons with small 
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virtualities to the amplitude decreases, and the suppression of double parton scattering is even 
stronger. Therefore the BFKL pomeron may be considered as a hard probe resolving the proton 
structure at virtuality \t\. 

The cancellation of the infrared divergencies is, in our approach, represented by introducing 
a physical cutoff so which is possible to interpret in terms of impact factors. Here we shall 
rather view it as a property of the confining QCD vacuum, in which propagation of gluons 
for distances larger than about 1 GeV^ 1 is suppressed |18|]. An approximate way to take this 



phenomenon into account is to modify the gluon propagator 1/k 2 — > l/(k 2 + so). 

In the high energy limit, the differential cross-section for quark-quark scattering is given by 



da qg (s,t) _ \A(s,t)\' 



dt 



16tt 



(1) 



where s is the squared invariant collision energy and t the momentum transfer. For the case of 
scattering by colour singlet exchange at large momentum transfer t, but s> \t\, the dominating, 



imaginary part of the amplitude is |L9| 



lmA(s,t = -Q 2 ) = J 



d 2 k 



% b {k,Q)^ ab {x,k,Q) 



(2tt)2 [(Q/2 + kf + s } KQ/2 - fc) 2 + s ] ' 



(2) 



Here, x = \t\/s, and k is defined so that the transverse momenta of the exchanged gluons are 
ki.2 = Q/2±k, where Q is the total transverse component of the momentum transfer t = —Q 2 . 
The impact factor for quark-quark scattering $Q 6 (fc, Q) = (5 ab /2N c )§ (k,Q) describes the 
coupling of the quark to two gluons with adjoint colour indices a, b in the colour singlet state, 
carrying the transverse momentum Q. The formulae for other relevant cases of quark-gluon and 
gluon-gluon scattering may be obtained from those given above by an appropriate modification 
of the colour charges. Beyond the LL approximation when the running of the QCD coupling g s 
(a s = <7 2 /47r) has to be taken into account, one has $o(fc, Q) = g s ((Q/2 + A;) 2 + s )g s ((Q/2 — 
k) 2 + So). In the LL approximation, however, a fixed value of g s will be used. $ ab is the result 
of the BFKL x-evolution from and is decomposed as § ab (x,k,Q) = (5 ab /2N c )$(x, k, Q). 
The function <3>(x, k, Q) satisfies the following BFKL equation 



$(x,fc,Q) = $o(&,Q) + 



3a s (fi 2 
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(3) 



This equation includes a treatment of the large next-to-leading corrections to the BFKL 
equation. They have a complicated structure [10], but their dominant part may be approxi- 



mately resummed to all orders by restricting the integration region in the real emission term 
f20f which is equivalent to resumming collinear divergencies in the NLL BFKL kernel 
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Generalising to the case of a non-forward configuration with Q 2 > the relevant limitation is 
H 

k' 2 <(k 2 + Q 2 /A)^ } (4) 

which follows from the requirement that the virtuality of the gluons exchanged along the chain 
be dominated by the transverse momentum squared. The constraint @ can be shown to 
exhaust about 70% of the next-to-leading corrections to the QCD pomeron intercept JTO, |20 



This constraint is embodied in eq. (|3]) by the step function Q((k 2 + Q 2 /4)x'/x — k' 2 ) which 
multiplies the real emission term. Another part of the non- leading corrections may be accounted 
for by using the running QCD coupling a s (fi 2 ) within the ladder, with the scale /z 2 related to 
the gluon virtualities. We use the one-loop QCD coupling a s with number of flavours Nf = 4 
and Aqcd = 0.20 GeV, with the scale fi 2 = k 2 + Q 2 /4 + s - Following our interpretation of the 
parameter sq, we have included it in eq. ([§) in a way which preserves gauge invariance. It is 
varied within the ranee 0.5 GeV 2 < s < 2 GeV 2 . Eq. (§ reduces to the standard LL BFKL 
equation after taking a fixed a s (we set ot s = 0.17) and substituting 0(...) — * 1. The equation 
is solved numerically in both cases; a summary of the numerical method and the adopted 



approximations is given in [21 



It is instructive to compare our approach with that of Mueller and Tang ||, who employ 
the exact solution of the non- forward LL BFKL equation |19) , given in terms of an infinite sum 
of contributions labelled by the conformal spin n. They retain only the n = eigenfunction of 
the BFKL kernel, which is the dominating contribution at high rapidity y and which has the 
anomalous dimension 7 = 1/2. Therefore their amplitude for elastic parton-parton scattering is 
convergent in the infrared and not sensitive to the details of the structure of the remnant. The 
remnant is only necessary to ensure the colour neutrality of the scattering hadrons. At small 
rapidities y, however, the full set of BFKL eigenfunctions has to be taken into account because 
of the completeness relation and the non-suppression of components with higher conformal 
spins at y — 0. 

The contribution of non-zero n at y — may be estimated by comparing the amplitude for 
two-gluon exchange 

including the infrared cutoff, with the Mueller- Tang amplitude 

19n 2 r v 2 

where u(v) = ^[2^(1) - ip(l/2 + iv) -^{1/2 - iv)\ and ip(z) = T'(z)/T(z). The ratio 

A 2g (t) ^ 1 



A MT (y = 0,t) 2 



log(l*l/*o) (7) 



is a large number for sufficiently large \t\. This ratio illustrates a mathematical property of the 
solutions to the BFKL equation, although at y — the high energy limit is strictly speaking a 
poor approximation. 
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The large logarithm in (|5|) is generated by asymmetric gluon configurations in which one of 
the exchanged gluons is soft with a small virtuality k\ ~ s > an d the other is hard with virtuality 
k 2 ~ Q 2 . The contribution of such two-gluon configurations to the colour singlet channel is 
suppressed for increasing y, because of an increasing probability of radiation from the hard 
gluon, which cannot be screened by the soft partner. This phenomenon is embedded into the 
virtual corrections to the BFKL kernel. The virtual corrections to the gluon propagation may 
be exponentiated to give a gluon reggeisation suppression factor {k\/Q 2 )^ asy ^ 2 ^ |T7|, After 
including this factor into the integral in eq. (|5|) we obtain an estimate of the non-asymptotic 
contribution of higher conformal spins 



. , . 16na 2 2tt 



9\t\ 3a s y 



-(3a a y/2n)' 



1 



e 2 s 



(8) 



where the e 2 factor was introduced in order to ensure A na (y — 0, t) + Amt(v = 0,t) = A^gif). 

A recent analysis of the contributions from higher conformal spins |22j shows that, in- 
deed, at smaller values of y they become dominant and account for the gluon reggeisation 
phenomenon. Thus, at moderate rapidities, such as what is now experimentally accessible, the 
original Mueller- Tang result is not sufficient for a proper description of the parton scattering 
amplitude at high \t\. Let us mention, that the presence of the scale Sq is reflected in the 
conformal representation by the appearance of an upper limit in the sum over conformal spins 



(i*i/ So ) i/4 m 



The ratio of the partonic differential cross-sections for colour singlet exchange and for leading 
order QCD quark-quark scattering by one-gluon exchange gives a naive estimate of the gap frac- 
tions before soft interactions and hadronisation are included. In fig. 0a we compare the numer- 
ical results of the LL BFKL equation, the Mueller- Tang equation (|6|), the non-asymptotic con- 
tribution corresponding to A na (y, t) as well as a simple sum of amplitudes A na (y, t) +A MT (y, t). 
The last guess is surprisingly close to the full solution. The non-asymptotic contributions dom- 
inate up to y = 6 despite the fact that the BFKL intercept 1 + ujq takes the rather high value 
of 1.45 which gives a very rapid increase of the Mueller- Tang n = component. 

Fig. ^|bc shows the corresponding naive gap fractions for quark-quark scattering resulting 
from the BFKL equation with non-leading corrections (|3D with sq = 1 GeV 2 and a running 
CK a (|t|) in dcr one _ gluon /dt. The decrease of the gap fractions continuing up to values of y larger 
than in the LL case is due to two effects. First, the non- leading corrections reduce the intercept 
and suppress the asymptotic, rising part. Second, the running coupling additionally enhances 
the contribution of asymmetric configurations (with one of the gluons soft) and presumably 
slows down the evolution of the reggeisation suppression factor with y. Furthermore, the gap 
fraction rises with E T ~ Q ~ ^J\t\ in spite of the fact that the running coupling was used. This 
is a specific property of the non-asymptotic part of the amplitude which still dominates even for 
y ~ 10. One should note that these dependences on y and Et of the gap fraction are in contrast 
to the predictions of the Mueller- Tang formula || |3|, |16|. This shows the importance of both 
subleading log(l/x) corrections and the non-asymptotic effects in y. Finally, it was checked 
that a variation of Sq between 0.5 and 2 GeV 2 does not influence significantly the shapes in 
and y. However, it does have a large effect on normalisation which increases (decreases) by a 
factor of about 1.5 when sq decreases from 1 to 0.5 GeV 2 (increases from 1 to 2 GeV 2 ) in the 
relevant range of y and Et- 
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3 Soft effects and Monte Carlo implementation 



In order to compare the perturbative results above with data, the parton level cross-section 
must be convoluted with the QCD evolved parton distributions fi/ p (x, \t\) and fj/p{x, \t\) of the 
incoming protons. In the absence of soft rescattering the cross-section for production of two 
jets initiated by partons i and j with a rapidity gap in between is given by 

where s is the invariant pp collision energy squared. The initial partons of species i,j carry 
momentum fractions x\^ (cf. fig. |l]a) and have colour charges Cjj = Cp for quarks and Ca for 
gluons. 

To consider rapidity gaps, which are defined as a region without final state particles, one 
must also consider all additional activity in an event. This means both pQCD processes, 
such as higher order parton emissions and multiple parton scattering, as well as soft processes 
including hadronisation. A simple way of doing this is to multiply the parton level cross-section 
(|9|) with a gap survival probability factor S 2 and sum over all parton species % and j. This 
probability may, however, depend on the subprocess with its associated colour field topology, 
or on the kinematics in the event. For example, in the case of gluon-gluon scattering or a 
larger momentum transfer, more partons are radiated and a smaller gap survival probability is 
expected. The best way to handle these complex processes is via Monte Carlo event simulation. 

We have therefore implemented the obtained cross-section for elastic parton-parton scatter- 
ing via colour singlet exchange, as well as the Mueller- Tang cross-section obtained from eq. @, 
as new hard subprocesses in the event generator Pythia. Higher order parton emissions from 
the incoming and outgoing partons are thereby included through conventional DGLAP parton 
showers and a model for multiple interactions (MI) |23j in terms of additional perturbative 



2^2 parton scatterings contribute to the underlying event complexity. Such a description 
is, however, incomplete since it only accounts for multiple exchanges of semihard gluons with 
Pt ~ Pt — ^ GeV 2 . Softer, non-perturbative interactions are modelled to take into account 
the remnants of the interacting hadrons, which are extended objects, and the hadronisation 



process through the Lund string model 



A successful approach to describe additional soft phenomena is given by the soft colour 
interaction (SCI) model @. Partons and remnants emerging from the hard process can exchange 
colour-anticolour, representing soft gluons with negligible momentum transfer. This leads to 
a modified colour string topology, producing a different final state after hadronisation. For 
example, a rapidity region without strings may arise resulting in rapidity gap events, and the 
SCI model does indeed reproduce the diffractive hard scattering processed observed at HERA 
25 and the Tevatron [pq| . Moreover, a colour octet cc quark pair may be turned into a colour 



singlet reproducing charmonium production [£7], [28[ . 

The soft colour interaction model is included in the Monte Carlo event generator together 
with the multiple interaction model in a way which is consistent with our previous study of 



diffractive hard scattering [ 26|| . As has been previously noted, SCI alone is not able to form 
sufficiently many gaps between jets when imposed on events with single hard gluon exchange 
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0. The reason is that the soft interactions cannot screen emissions of gluons in the hard 
subprocess. Furthermore, contrary to the perturbative semihard gluons in eq. the SCI 
contribution is not enhanced by a loop integration. In the present case of gaps between jets, 
the SCI has the effect of rearranging strings to span across the parton level gap from the BFKL 
singlet exchange and therefore it reduces the gap survival probability. 



4 Comparison with data 

Gaps between jets have been observed in pp collisions at y/s = 1800 GeV by the D0 || and 
CDF 0] collaborations. The results are presented in terms of the fraction of two-jet events that 
have a rapidity region \rj\ < 1 with no particles in between the jets. The jets are reconstructed 
using a cone algorithm and have transverse energy defining different samples. D0 has a 
Iow-Et sample with 15 < Et2 < 25 GeV, and a high-.&r sample with Et2 > 30 GeV, where 
Et2 is the second highest jet transverse energy. CDF has one sample with Et > 20 for both 
jets. 

The same criteria and analyses as for the data can be applied to the Monte Carlo events 
generated with the parton-level BFKL cross-section (including non-leading corrections). The 
gap fraction is given by the ratio of hard colour singlet exchange events to the standard LO 
QCD processes. The effects of the underlying event leading to the destruction of rapidity 
gaps are studied in detail and illustrated by three different treatments. The first and simplest 
method is to multiply the generated cross-section for gaps at the parton level with a constant 
survival probability S 2 taken as a free parameter. The second method is to use the MI model 
(in its most sophisticated version given by setting MSTP(82)=4 in Pythia) followed by standard 
hadronisation. The third and most elaborate method is to use both MI and SCI before applying 
string hadronisation. The value of the infrared regulator of the MI model pt has here been 



increased to avoid double counting of soft interactions, as described in |26 . 

In fig. |3] we compare the model results with the D0 data. First, one can notice that the 
Mueller- Tang approximation with constant gap survival probability gives a wrong Ej- depen- 
dence of the gap fraction as well as a somewhat too large slope of the rapidity dependence. 
Going beyond the large- y approximation and including the non-leading BFKL corrections, gives 
an improved result that can reproduce the data quite well. Without our most elaborate model 
of the underlying event, however, one needs an overall renormalisation in terms of a gap sur- 
vival probability. Since no such factor is needed when both the multiple interaction model and 
the soft colour interaction model are included, it seems that these models account properly 
for the gap survival probability. The normalisation is here rather sensitive to the value of the 
probability parameter P for a soft colour exchange between partons. For example, the gap 
cross-section decreases by a factor two when increasing P from 0.20 to 0.30. In contrast, the 
shapes of the y and ^-distributions are not sensitive to variations of P. We have used the 
value P = 0.25, which is within the range allowed by the previous analyses of the SCI model 
discussed above. 

One should note that the solution to the BFKL equation with non-leading corrections gives 
a shape of the Et distribution (fig. ^a), which is stable against variations of the model of the 
rescattering process in the underlying event. This gives strong support for the BFKL approach 
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and shows that the Mueller- Tang approximation is not valid in the kinematical range of the 
Tevatron. The Arj dependence is much more sensitive to the treatment of rescattering processes. 
The inclusion of multiple interactions slightly increases the slope of the rapidity dependence 
and an even stronger increase is observed when the SCI model is also turned on. 

The CDF data for the gap fraction as a function of Arj are shown in fig. |j. A discrepancy 
with the BFKL+MI+SCI model can be observed in the region Arj > 5.5. The CDF data 
show a different tendency than the D0 data, but because of the large errors it cannot be 
claimed that the data are incompatible. The BFKL result with constant gap survival probability 
S 2 = 3% falls between these data sets and therefore gives a reasonable overall description. 
Recall however, that in this case S 2 is arbitrary. 

In fig. [|b, this model is decomposed into the contributions from the different hard scattering 
subprocesses. Quark-gluon scattering dominates except at the largest Arj where quark-quark 
scattering contributes equally, whereas gluon-gluon scattering is only important at the lowest 
A77. These relative contributions can be understood, since both gluon and quark initiated 
processes contribute to the cross-section (|9|) in these jet samples where the initial partons 
typically have a fraction x ~ 0.1 of the beam momentum. The gluon contribution in the 
HCS exchange is enhanced in relation to the quark one by (Ca/Cf) 2 , the ratio of the colour 
charges squared. This partially compensates for the smaller value of the gluon distribution in 
the probed region of x. For At] ~ 4, x may be smaller and one has more gluon scattering, 
whereas for Arj ~ 6 quark scattering dominates. 

The increase of the rapidity slope of the gap fraction caused by the SCI model can now 
be traced. The colour string topology in the Lund model depends on which parton that was 
scattered. A scattered valence quark leaves a diquark remnant as endpoint of a string, whereas 
a scattered gluon leaves a diquark plus a quark stretching two strings to the remnant. Thus, 
with a constant probability for colour exchange within parton pairs in the SCI model, there is 
a larger probability for string reconnections with the more partons and strings available in case 
of gluon scattering. For the hard colour singlet exchange this means a larger probability to get 
a string across the parton level gap in case of gluon scattering. The gap survival probability is 
therefore smaller for gluon-gluon scattering, being more frequent at small Arj, and largest for 
quark-quark scattering at larger Arj. 

Finally, let us recall that the gap fraction is defined relative to the total rate of high- 
er jet pairs with the same rapidity separation. This corresponds to production of Mueller- 
Navelet jets ||29|| , which is not well understood when the same Ex-cut is imposed on both 
jets. The experimental results suggest a rather steep dependence of the Mueller- Navelet cross- 
section on the rapidity separation between the jets [|30|. Recent theoretical considerations [31 



indicate a BFKL-driven increase of the cross-section for production of Mueller-Navelet jets 
(with Et cut — 20 GeV) by about 50% when Arj increases from 4 to 6. However, the equal cut 
on Et makes this observable infrared unsafe, i.e. , unstable against variations of the treatment 
of higher order QCD corrections. It is still fair to say that both theory and experiment give an 
expectation of an increase with rapidity of the denominator of the gap fraction. This would, 
of course, have some impact on the Arj shape of our predictions. We have used Pythia with 
conventional 2 — > 2 QCD matrix elements for the denominator in the gap fraction. Inclusion 
of the suggested increasing tendency would affect our model results and reduce their A^-slope. 
This would bring the results obtained with MI+SCI in better agreement with the CDF data in 
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fig. ^a, whereas the simpler model with constant gap survival probability would get a negative 
slope that spoils its agreement with data. Lacking a reliable theoretical understanding of the 
Mueller-Navelet jet production mechanism, we leave this problem open for future studies. 

5 Conclusions 

We have investigated the production of high-.&r jets separated by a central rapidity gap as 
observed in pp collisions at the Tevatron. The data can be understood in terms of a model 
based on hard parton scattering through the exchange of a colour singlet gluon ladder described 
by the BFKL equation. The two-gluon exchange amplitude with resummed virtual corrections 
plays a crucial role in the Tevatron kinematic range, where the energy is not asymptotically 
large. This makes the previously used Mueller- Tang approximation inadequate. The account 
for non-asymptotic effects in y and non-leading corrections in the BFKL formalism improves 
the theoretical understanding and results in distributions in Ej- and Arj of the parton level 
gap fraction in much better qualitative agreement with the data. Our implementation of 
the partonic BFKL process in the Pythia Monte Carlo provides a detailed treatment of soft 
dynamics giving a gap survival probability that varies event by event. The full event simulation 
is particularly important for the rapidity gap observable defined by the final state properties. 

A detailed comparison with the D0 and CDF data shows an overall good agreement. In 
particular, the model is able to reproduce the magnitude of the gap fraction. The E T distribu- 
tion of the gap fraction coincides with the experimental results. This conclusion holds both at 
the partonic and hadronic level regardless of the treatment of soft effects. On the other hand, 
our result for the rapidity dependence of the gap fraction is somewhat uncertain. This problem 
is not related to the hard colour singlet exchange but to an unresolved issue of rapidity depen- 
dence of the cross-section for the production of Mueller-Navelet jets. Within this uncertainty, 
a reasonable description is obtained for the Arj distribution of the gap fraction. 

In conclusion, the BFKL approach to hard colour singlet scattering gives a good under- 
standing of the jet-gap-jet events observed at the Tevatron. This demonstrates the importance 
of BFKL dynamics in QCD. 
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Figure 2: Naive parton level gap fractions [da S i ng i et (y, t) / dt]/[da one - g \ non (y, t) / dt] in quark-quark 
scattering via colour singlet exchange and via leading order one-gluon exchange, versus the 
rapidity difference y of the two scattered partons and their transverse energy Et- (a) Leading 
order BFKL solution (fixed a s = 0.17) for Et=25 GeV compared with the Mueller- Tang 
asymptotic approximation to it, eq. ([]), the reggeised gluon approximation eq. (^), and the 
sum of the two. (b,c) BFKL solution with non-leading corrections. 
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Figure 3: Fraction of jet events having a rapidity gap in \r)\ < 1 between the jets, versus (a) the 
second highest }et-E T and (b,c) the rapidity separation Ar/ between the jets for low and high 
jet-E T (15 < E T2 < 25 GeV and E T2 > 30 GeV). D0 data compared to the colour singlet 
exchange mechanism based on BFKL equation with non-leading corrections with the underlying 
event treated in three ways: simple 3% gap survival probability, multiple interactions (MI) and 
hadronisation requiring a 15% gap survival probability, MI plus soft colour interactions (SCI) 
and hadronisation with no need for an overall renormalisation factor. Also shown is the Mueller- 
Tang (MT) calculation with an 11% gap survival probability. 
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Figure 4: Gap fraction as in fig. |3], but normalised to the average gap fraction, versus the rapidity 
separation A?7 between the jets, (a) The BFKL colour singlet exchange model with (solid 
curve) and without (dashed curve) underlying event simulation, calculated with conditions as 
for the shown CDF data ||. Rescaled D0 data from fig. |3| are included for comparison, (b) 
BFKL result for the contributing subprocesses of qq, qg and gg scattering, but without multiple 
interactions and soft colour interactions. 
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